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5^ ■ Abstract. We present the numerical results for low temperature behavior of the transverse- 

T^ ' field Ising model on a frustrated checkerboard lattice, with focus on the effect of both quantum 

and thermal fluctuations. Applying the recently-developed continuous-time quantum Monte 
C^ ' Carlo algorithm, we compute the magnetization and susceptibility down to extremely low 

temperatures while changing the magnitude of both transverse and longitudinal magnetic fields. 

Several characteristic behaviors are observed, which were not inferred from the previously- 
' O ' studied quantum order from disorder at zero temperature, such as a horizontal-type stripe 

C I ordering at a substantial longitudinal field and a persistent critical behavior down to low 

>-» . temperature in a weak longitudinal field region. 



>, 
^5 ' !• Introduction 

Frustration due to the geometrical structure of systems has been one of the major topics 

in condensed matter physics. The interest covers a wide range of fields, such as frustrated 

r~^ . magnetism [T], proton ordering in hydrogen-bonded systems [21 [3], and anomalous transport 

^^ I and multiferroics in transition metal oxides [H El [B] . The geometrical frustration often prevents 

the system from selecting a unique ground state, instead gives rise to classically degenerate 

ground states; i.e., a macroscopic number of different configurations of the system variables lead 

to the same ground state energy at the classical level. As a result, the ground state remains 

S^ . to be disordered and bears finite residual entropy. Such macroscopic degeneracy provides a 

H I fertile ground for various peculiar phenomena as it makes the system to be extremely sensitive 

to perturbations, such as remnant interactions, external fields, and fluctuations. 

Among such degeneracy-lifting mechanisms, quantum and thermal fluctuations have attracted 
particular interest. These fluctuations sometimes lift the ground-state degeneracy and induce 
some particular ordering, as known by order from disorder [7j. There, the ordering is selected 
from the manifold to maximize the associated entropy (zero-point entropy in the case of quantum 
fluctuation) . This entropic effect plays a decisive role at low temperatures (T) and causes many 
fascinating phenomena in the frustrated systems. 

One of the minimal models for studying the order-from-disorder phenomenon is the 
transverse-fleld Ising model (TIM). In the absence of geometrical frustration, in general, the 
model develops a long-range order at low T, while the ordering is suppressed and a quantum 
paramagnetic state is induced by the transverse field. When the classical ground state is 



macroscopically degenerate under strong frustration, quantum fluctuations by the transverse 
field and/or thermal fluctuations by temperature can induce a particular ordering via the order- 
from-disorder mechanism. Efl'ect of thermal fluctuations has been studied mainly in the absence 
of the transverse field, i.e., for pure Ising models without the transverse field [Il[8]. Meanwhile, 
effect of quantum fluctuations has also been studied. For example, it was shown that TIM on a 
variety of frustrated lattices exhibit several nontrivial behaviors at T = 0, such as a bond ordering 
and Kosterlitz-Thouless transition OITO]- In general, the thermal and quantum fluctuations do 
not necessarily lead to the same effect, and the relation between them is of particular interest 
to explore yet another order-from-disorder phenomenon. 

In this contribution, we present our numerical results for the order-from-disorder phenomena 
in TIM in a wide range of temperature and the transverse/longitudinal magnetic fields. We 
consider the TIM on a two-dimensional checkerboard lattice and study its low-T physics by 
a sophisticated quantum Monte Carlo (QMC) technique. As a result, in the intermediate 
longitudinal-fleld regime, we find instability toward a Neel order, in accordance with the previous 
report [10]. On the other hand, for both the weaker and stronger longitudinal fields, our results 
indicate different behaviors from the previous report. Under the weak field, featureless magnetic 
susceptibility is observed down to extremely low T, which implies a very weak proximity effect 
to the predicted Neel ground state or a possibility of another state in the T = limit. On the 
contrary, in the strong field, the system indicates instability toward an unexpected horizontal- 
type stripe ordering. 

2. Model and method 

To investigate the effect of quantum and thermal fluctuations on the classically degenerate 
ground state in frustrated systems, here we focus on the frustrated checkerboard Ising model 
with transverse and longitudinal magnetic fields. The Hamiltonian is given as: 

H = jY,sls^ + rJ2sf + h^si (1) 

where s" is the Pauli spin operator at site i, and J is the Ising interaction between nearest- 
neighbor sites {i,j) on the checkerboard lattice (see Fig. [1]); F and h are the transverse and 
longitudinal fields, respectively. In the following, we consider the antiferromagnetic case, J > 0, 
and set the energy scale by J = 1. 

When F = /i = 0, the model remains to be disordered down to zero T [8]. The lowest energy 
is achieved by enforcing a simple local constraint — a two-up two-down local spin configuration 
for all the plaquettes with crisscrossing interaction, similar to the so-called ice rule in water ice 
[2]. Consequently, a macroscopic number of different spin configurations give the same lowest 
energy: The ground state is macroscopically degenerate and the residual entropy is estimated to 
be ~ Y log I (A'^ is the number of spins). The spin correlation, however, is critical in the sense 
that it decays algebraically as a function of distance. The situation is unaltered as h increases 
up to /i = 2J at F = 0. (For larger h, the ground state consists of "three-up one-down" 
configurations, which is also macroscopically degenerate.) When the quantum fluctuation sets 
in by switching on F, it was predicted, from perturbative considerations in T/h, that a Neel 
order is induced in the < h < 2J region at T = [TU] [see Fig. [U^a)]. Hereafter, we focus on 
the 0<h<2J region. 

To investigate the low-T behaviors of this model with taking account of both quantum 
and thermal fluctuations, we conducted numerical calculations by QMC method based on the 
path-integral approach [IH 112] . Here, we employed the recently-developed continuous-time 
algorithm [13] to approach extremely low T. At low T, the MC sampling suffers from the slow 
relaxation due to the severe frustration. To accelerate the relaxation process, we applied the 




Figure 1. Schematic pictures of 
tlie transverse-field Ising model 
on tlie clieckerboard lattice, (a) 
Neel type order and (b) horizon- 
tal stripe order. 

loop-flip algorithm developed for frustrated Ising-type spin systems |14U15j . In our calculations, 
the loops are formed in the real space at a particular imaginary-time slice, and flipped together 
with all the spin variables along the imaginary-time direction, i.e., on the torus defined by the 
loop in space and time. In addition, the replica exchange method was used to further suppress 
the slowing down [16]. For the present system, the replicas are exchanged along a constant-F/T 
line because the Boltzmann weight depends on the number of domain walls in the imaginary-time 
direction which is proportional to F/T. 

Calculations were done for N = Ax L"^ site systems with L = 12 under the periodic boundary 
conditions. Observables were computed typically for 120000 samplings after 30000 steps of 
thermalization. Fhe results were divided into eight bins to evaluate the statistical errors. 

3. Results and discussions 

Figure [2] shows the QMC results for the susceptibility at a relatively weak longitudinal field 
h = 0.4 with varying F. Figure [2]||a) presents the susceptibility for the Neel- type order parameter 

[Fig. (Ha)], 

XNccl = ^(("^Nccl) - (kl^Nccll)^) ; mNeel = - J]] <(-!)'" ^'^ (2) 

i 

while Fig. El^b) is that for the horizontal-type stripe order parameter [Fig. (D^b)], which we 
consider as another candidate of order from disorder, 

Xstripe = ^(('^stripc) " Ktripc)') ; ^^^ripe = {^Yl ^'("l)'') + (]v ^ 4(-l)'") • (3) 

i i 

Note that this stripe state belongs to the two-up two-down manifold and is different from the 
diagonal one discussed in the three-up one-down manifold for 2J < h < 6J [10]. Figure [2]|^c) 
shows the local correlation parameter given by /? = jf'^pfip), where the sum runs over all 
plaquettes. f{p) is a function defined for each plaquette p as f{p) = 1 for "^i^pSf = and 
otherwise f{p) = —3/5; p — )• 1 when all the plaquettes satisfy the two-up two-down spin 
configuration, and p ^>- when the spins are completely disordered, i.e., as T — )• oo and/or 
F — )■ oo. As shown in Figs. [2]^a) and[2]J^b), for weak F < 0.5, both the susceptibilities increase 
with decreasing T in Curie-law like behavior, i.e., x oc T~^. A shoulder-like feature observed 
at r ~ 0.3 in Xstripc corresponds to a formation of the ice-rule type local configurations, as 
indicated by the saturation of p in Fig. EJJ^c). These behaviors suggest that the system is in 
the correlated regime with satisfying the local constraint, but remains critical, without choosing 
either Neel or stripe ordering, presumably because of the strong frustration. We also calculate 
the momentum dependence of the susceptibility, and confirm that it remains featureless without 
showing any peak in the momentum space. It is surprising that the critical behavior is robustly 
observed down to very low T ~ 0.01. This implies that the order-from-disorder mechanism is 
extremely weak or ineffective in this region. On the other hand, the results for F > 0.6 show 
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Figure 2. T dependence of (a) XNeeb (b) Xstripe, and (c) local correlation parameter 
h = 0.4 for the transverse field ranging from T = 0.11 to 0.61. 
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saturating behavior at low T. In this region, T is strong enough to disturb the ice-rule type 
configurations [Fig. [2jc)], and the ground state turns into a quantum paramagnetic state. 

With increasing the longitudinal field, a contrasting behavior between the two susceptibilities 
shows up in the vicinity of phase transition to the quantum paramagnetic state. Figure [3] shows 
the results at h = 0.9. In the region F < 0.3, both xnccI and Xstripe diverge toward T = 0, in 
a similar manner to the results for F < 0.5 in Fig. [2j On the other hand, before entering the 
quantum paramagnetic regime for F > 0.5 where both xnccI and Xstripe saturate, a qualitatively 
different behavior is observed: At F ~ 0.40, Xstripc begins to deviate from the Curie-law like 
behavior and tends to saturate, while XNeei remains to diverge (indicated by the arrows in Fig. [3]). 
Fhis implies that the stripe-type fluctuation is suppressed and the Neel-type ordering is likely 
favored in the low-T limit. Similar behavior is observed in a finite region on the verge of the 
quantum paramagnetic phase for 0.5 < h < 1.0 (see the phase diagram in Fig. [5] below). 
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Figure 3. T dependence of (a) XNccIi (b) Xstripci and (c) local correlation parameter p at 
h = 0.9 for the transverse field ranging from F = 0.11 to 0.61. Fhe arrows in (a) and (b) 
indicate a contrasting behavior between XNeei and Xstripc in the intermediate F region. 



For larger longitudinal fields, yet another behavior appears. For h > 1.0, the MC relaxation 
becomes very slow at low T even with the use of the loop-flip algorithm and the replica-exchange 
technique. The situation is demonstrated in Fig. [H The results show the relaxation of the order 
parameters rriNeei and TTistripe at /i = 1.6 along a constant T/T = tan(85°) axis; the data are 
measured for 50000 MC steps after particular thermalization steps of A'therm starting from the 



initial configuration witli corresponding perfect order: Figures Sl^a) and Hl^b) show tlie results 
starting from the perfectly-Neel-ordered state and the perfectly-horizontal-stripe-ordered state, 
respectively. As shown in Fig. HI JTiNeci ~ after A'therm = 10^ thermalization independent of the 
initial state. On the other hand, the w-stripe shows strong dependence on the initial state even 
after A'^thorm ~ 10^ steps. A possible origin of the freezing behavior is strong first order transition 
to the horizontal-stripe-ordered state. Although the extremely slow relaxation prevents us from 
obtaining numerically-converged results, the results suggest that the system has an instability 
toward the horizontal stripe ordering at low T in this high-/i region. 
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Figure 4. T dependence of ttinccI and rristripc at /i = 1.6 and T/T = tan(85°) for 
varying thermalization steps A^therm- (a) is the results starting from the perfectly-Neel-ordered 
configuration and (b) is from the perfectly-horizontal-stripe-ordered one. 



Summarizing the results above, we deduce a 'phase diagram' in Fig.[Sl which indicates the low- 
T instability anticipated from the MC data. In the strong longitudinal field region, 1.0 ^ h < 1.8 
and r < 0.2, our MC results indicate an instability toward the horizontal-stripe-type ordering 
below some particular temperature, as exemplified in Fig. Sj^b). Although we could not directly 
confirm the phase transition because of the extremely slow MC relaxation, we expect this phase 
to appear at a finite temperature. In the intermediate longitudinal field region, 0.5 ^ h ^ 1.0 
and 0.3 < F < 0.5, the results appear to favor the Neel-type instability. The Neel temperature 
appears to be extremely low, less than T = 0.01. In the remaining low-h region for F < 0.5, 
we could not observe any sign of clear instability toward either Neel or stripe type, at least for 
T > 0.01; the system is strongly correlated but remains to be critical under the frustration. 

Let us compare the result to the T = argument by Moessner and Sondhi |10] . In the previous 
study, by a perturbation in T/h, it was predicted that a Neel-type order occurs in the entire 
region of our consideration. Our result in Fig. [5] shows a similar tendency in the intermediate 
h, where the Neel-type order appears to be favored. However, we could not detect any sign of 
a particular ground state in the lower-/i region, although the Neel-type order was expected to 
extend down to /i — > in the schematic T = phase diagram in the previous study [10]. This 
apparent disagreement can be due to very weak degeneracy lifting by the order from disorder 
for the expected Neel ground state. Another possibility is that, for large T/h, the perturbative 
argument in T/h no longer holds and a different phase emerges in the region. In order to identify 
the nature of the system in this region, further analysis down to lower T in larger system sizes is 
necessary. On the other hand, under relatively strong field 1-0 < /i < 1.8, we detected instability 
toward another phase from the slow relaxation behavior, namely, the horizontal-stripe phase. 
This might be related with the competition in the ground states between the two-up two-down 
Neel order and three-up one-down diagonal-stripe order at /i ~ 2. Our preliminary calculations 



show a finite-size effect in this regime; the horizontal-stripe type instabihty appears to extend to 
shghtly lower-/i region in larger size systems. More detailed results will be discussed elsewhere. 

Figure 5. 'Phase diagram' for the model in 
Eq. ([1]) as functions of transverse field T and 
longitudinal field h. The results indicate the low- 
T instabilities which are deduced from the QMC 
j!j 1 _ ! ! I _ results down to T ~ 0.01. The blue squares 

(orange circles) show the region in which the 
Neel(horizontal stripe)-type instability shows up. 
The green triangles show the region where both 
XNcel and Xstripe remain to show the Curie-law 
like divergence down to the lowest T. The area 
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 without symbols for T > indicates the quantum 
Y paramagnetic region. 



4. Summary 

To summarize, we have investigated the effect of quantum and thermal fluctuations on the 
thermodynamics of the frustrated checkerboard Ising model with transverse and longitudinal 
magnetic fields. We have examined the magnetic instabilities induced by order-from-disorder 
mechanism by employing a continuous-time quantum Monte Carlo method with the loop-flip 
update and replica exchange algorithm. We identified characteristic behaviors, which were not 
inferred from the previous study for the quantum order from disorder in the ground state. One 
is the instability toward a horizontal-stripe-type ordering in a relatively high longitudinal-field 
regime. This is distinguished from the diagonal stripe predicted for a higher field. The other is no 
clear indication of dominant ordering or fluctuation down to T ~ O.OIJ in a relatively low- field 
regime; the susceptibility shows a persistent Curie-law-like divergence at all the momenta. This 
could be due to the emergence of a novel state or simply by surprisingly weak order-from-disorder 
effect in the Neel ground state expected in the previous theory. 
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